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Abstract. We perform an analytical study and a Monte Carlo (MC) analysis of the main features for microlensing 
events in pixel lensing observations towards M31. Our main aim is to investigate the lens nature and location of 
the 14 candidate events found by the MEGA collaboration. Assuming a reference model for the mass distribution 
in M31 and the standard model for our galaxy, we estimate the MACHO-to-self lensing probability and the event 
time duration towards M31. Reproducing the MEGA observing conditions, as a result we get the MC event 
number density distribution as a function of the event full-width half-maximum duration t\n and the magnitude 
at maximum Rmax. For a MACHO mass of 0.5 Mq we find typical values of tj/ a ~ 20 day and i? max — 22, 
for both MACHO-lensing and self-lensing events occurring beyond about 10 arcminutes from the M31 center. 
A comparison of the observed features (ti/a and -Rmax) with our MC results shows that for a MACHO mass 
> 0.1 Mq the four innermost MEGA events are most likely self-lensing events, whereas the six outermost events 
must be genuine MACHO-lensing events. 
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1. Introduction 

Gravitational microlensing has become since about a 
decade a robust tool for analyzing the galactic structure 
and for gaining infor mation about the dark mass com- 
pone nt in our galaxy ijAlcock et all Il993t lAubourg et all 
1199.1) . Several hundreds of microlensing events have been 
detected so far towards the galact ic bulge, the spiral 
arms, and the M agellanic Clouds 1 Alcock et all [2000; 
I Jetzer et all 12004^ . 

Recently, pixel lensing observations towards M31 and 
even M87 have been undertaken and some microlens- 
ing events have been found. The MEGA collaboration 
ijde Jong et all l2004[) has reported the detection of 14 
candidate events towards M31 by using Isaac Newton 
Telescope (INT) data. Two of these events have been 
previously reported by the POIN T- AGAPE collaboration 
( Paulin-Hcnri ksson et al 1 1200.1) . which recently has pre- 
sent ed a high-threshold analys is of the full 3 years data 
set l|Calchi Novati et all . 120051. This analysis shows that 
the observed signal is much larger than expect from self- 
lensing alone and that some fraction of the halo mass must 
be in form of MACHOs. Other collaborations have also re- 
ported preliminary r e sults f or pix e l lensing toward s M31 
ilCalchi Novati et all I2002L 1200.1 llJglesich et all l2004t 



didate s. A m ap of the MEGA events is reported by de Jong 
et al. l|2004|) . A preliminary analysis by the MEGA col- 
laboration indicates that the events located in the outer 
part of M31 are consistent with being due to halo lens 
objects, whereas the innermost ones are most likely due 
to self-lensing. The aim of this paper is to perform a more 
complete analysis by using a Monte Carlo (MC) program 
with the purpose of using all available information by the 
MEGA collaboration, which are summarized in Tab.^ to 
characterize the nature of the lenses. 

The paper is organized as follows: in Section 2 we 
present the calculation of the microlensing rate towards 
M31 and we remind the pixel lensing basics. In Section 
3 we give some details on the adopted M31 and Milky 
Way (MW) mass distribution models as well as the stellar 
and mass functions. The analytic and MC results and a 
comparison between the two are given in Sections 4-6 and 
conclusions are presented in Section 7. 



2. Microlensing rate 

T he differential number of expected m icrolensing events 
is jPe Ruiula et allll99ltlGriestlll99l|) 



nuaicni iNovat i et ail izuua. i^uudt lug. 
iRiffeser et all I2001 Ijoshi et all 120041) . 

In this paper we consider in particular the MEGA re- 
sults which provide the largest sample of microlensing can- 



dN cv = N*t ohs dT, 



(1) 



where is the total number of monitored stars during 
the observation time i Q bs- The differential rate dT at which 
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Table 1. For the 14 MEGA events we give the position, 
the full- width half-maximum duration ti/2 and the magni- 
tude at maximum R max . The coordinate system we adopt 
has origin in the M31 center and the X axis is oriented 
along the M31 disk major axis (see also Fig. QJ. 



MEGA 


X 


Y 


i obs 
1/ z 


r> o 






arcmin 


arcmin 


day 


magn 


1 


-4.367 


-2.814 


4.20 ± 4.30 


22.2 


± 


1.1 


2 


-4.478 


-3.065 


4.60 ± 0.60 


21.6 


± 


0.3 


3 


-7.379 


-1.659 


2.60 ± 2.20 


21.8 


± 


1.2 


4 


-10.219 


3.420 


29.10 ± 1.00 


22.8 


± 


0.2 


5 


-19.989 


-13.955 


9.40 ± 4.10 


22.9 


± 


0.8 


6 


-21.564 


-13.169 


22.90 ± 0.70 


22.6 


± 


0.2 


7 


-21.163 


-6.230 


21.60 ± 0.70 


19.3 


± 


0.2 


8 


-21.650 


7.670 


27.40 ± 0.90 


22.7 


± 


0.2 


9 


-33.834 


-2.251 


3.80 ± 1.60 


21.8 


± 


0.8 


10 


-3.933 


-13.846 


46.80 ± 4.40 


22.2 


± 


0.3 


11 


19.192 


-11.833 


2.00 ± 0.30 


20.5 


± 


0.2 


12 


29.781 


-5.033 


131.00 ± 9.40 


23.2 


± 


0.3 


13 


22.072 


-22.022 


22.80 ± 3.80 


23.3 


± 


0.3 


14 


19.348 


-29.560 


28.10 ± 1.40 


22.5 


± 


0.2 



a single star is microlensed by a compact object is given 

by 



dT 



n\(x, p)dpf(vi±)d 2 v\±d 3 o 
Jt 



(2) 



where the numerator on the right hand side is the number 
of lenses with transverse velocity in d 2 v\± = v\± dv\±_ d/3 
around vi±, located in a volume element d 3 x = dx dy dz 
centered at the position x of the microlensing tube. Here 
v\± is the component (in the rest frame of the Galaxy) 
of the lens velocity orthogonal to the line of sight to the 
source star, n\(x,p) the lens number density (per unit of 
volume and mass) and f(v\±) the lens tranverse velocity 
distribution, p is the lens mass in solar units. 

In evaluating dT we must take into account that the 
source stars are not at rest but have a transverse ve- 
locity v s ± as well, which can be splitted into a random 
component and a component which describes the ordered 
rotation (if present) of the galactic component consid- 
ered, namely v s ± = i> s j_ ir an + ^s±,rot- Using a Maxwellian 
distribution with 2-dimensional velocity dispersion o~ s to 
describe the random velocity component, we get for the 
transverse velocity distribution g(v s ±) 



g(v s ±)v s xdv s ±d(p = — i-^e 



v s ±dv s xdcp, (3) 



where ip is the angle between v s ± and i> s ±,rot- 

We take also into account the transverse rotation ve- 
locity of the observer, v@ ± rot (here we consider the ob- 
server co-moving with the Sun) by considering that the 
microlensing tube is moving with a transverse velocity 



1 



£>os 



v @X,rot 



A 

Dr 



1 

-v s ±_ , 



(4) 



where D os and D a \ are the source and lens distances from 
the observer, respectively. 

As for the source stars we split the lens transverse ve- 
locity into a random component and an ordered rotation 
component. Moreover, by taking into account that the mi- 
crolensing tube moves with velocity i>t_L, it follows that 
V\J_ = ^l_L, ran + V\x,rot — Vtx- Accordingly, assuming for 
v\±.ran a Maxwellian distribution with 2-dimensional ve- 
locity dispersion er; , we obtain using also the definition for 
g{v s ±) given in eq.J2J) 



f(vi±)vi±dv\x d(i 



1 



{-K(J s <Ji) 



dtp x 



o 



2 

v s ±dv s ± e "i v l± dvi± d(3 , (5) 



where w± = v\j_ trot — v t ±_ and (3 is the angle between v\± 
and wj_. 

We can then write the volume element, d 3 x, as 



d 3 x = (vi± ■ h) dt dS = vi± cos 9 dt dl dD a \ 
= cos6> dt i?E du t h da dD Q \, 



(G) 



where 9 is the angle between v\±_ and the normal, n, to 
the lateral superficial element, dS — dldD Q \ 1 of the mi- 
crolensing tube, with dl = R^du^da being the cylindrical 
segment of the tube (wth is the threshold value for the im- 
pact parameter). Note that, if a is taken to be the angle 
between h and w± , it follows that 9 — a + (3, so that for a 
constant value of /?, da = d9. Therefore, the microlensing 
differential event rate becomes 

dT = ni(x,fj,)diJ 1 f(vi±)v\ ± dv\±df3cos9REdu t i 1 dadD \ . (7) 

We assume, as usual, that the mass distribution of 
the lenses is independent of their position in M31 or in 
the Galaxy (factorization hypothesis). So, the lens num- 
ber density n\(x, /i) can be written as ijjetzer et all 12002^1 



n\(x,ii) 



Po 



V'o(m) , 



(8) 



where po is the local mass density in the Galaxy or the cen- 
tral density in M31, ipo([i) the corresponding lens number 
density per unit of mass and the normalization is 



/'up 



V>o(m) M dp = 

Mq 



(9) 



Here p m i n and p up are the lower and the upper limits 
for the lens masses (see Subsect. 3.3). Accordingly, the 
microlensing event rate is given by: 



rpos) 



y — J dp p 1/z ip (n) J duth x 



o V 

2lT /"OO /> + 7r/2 

d/3 / d«ix/(«ix)vu / 

n JO J-tv/2 



Po 

cos 9d9 , (10) 
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Uthreshold map 




.13 

.14 



Fig. 1. The (ut(x, y))^ RG contour plot towards the M31 
galaxy is shown with the position in the sky plane of the 
14 MEGA candidate microlensing events. From the inner 
to the outer part of the figure (ut(x, y))^> RG increases with 
lines referring to the values 0.02, 0.04, 0.06, 0.08 and 0.10, 
respectively. Outside the 0.10 contour line (uT(x,y)),p RG 
is almost constant. Note that the positions of the events 
1 and 2 are almost on top of each other. The same is true 
as well for all other figures. 



where the integration on 9 is performed between — 7t/2 and 
+tt/2, since only lenses entering the microlensing tube are 
considered. After integrations on 9 and f3 we finally get 



r(D os ) = 2a\UT 



4GMr. 



' ' ' dn /i 1/2 Mp) x 



D,, 



dD n 



1 D ol (D os - D ol ) Pl (D ol ) 



Po Jo 



P{z)dz , (11) 



where z = V\x/o\ and the function P(z) is given by 



P{z) 



2e~ a 



271 



dip / ye~ {y -iaycos V + v ) dy x 

i) Jo 

z 2 e~ z2 1 (2r]z) . (12) 



Here a(D os ) = v s j_ <rot /cr s , y = v s± /a s , r){D OB ,y,ip) = 
w±/ai and Io{2rjz) is the zero-order modified Bessel func- 
tion of the argument 2r\z. 

Note that x( z ) — z^ 1 P(z) is the dimensionless 
form of cq. © and that it is properly normalized as 
it can be easily verified by using twice the relation 
J °° dt exp{-t 2 ) t h{2qt) = exp{q 2 )/2. 

To take into account the source distribution in the M31 
bulge and disk, eq. (|llfl has to be integrated not only over 
the distance of the lenses but also over the distance of the 
sources. Accordingly, the microlensing rate becomes 



/ °° Ps (D 0S )T{D 0S )dD 
J Q °° p s (D os )dD os 



(13) 



where x and y are coordinates in the plane orthogonal to 
line of sight, p s is the source mass density (which is the 
sum of the sources in the M31 bulge and disk). 

Moreover, we compute the average Einstein time 
(which depends on the line of sight position given by the 
coordinates x and y) as 

la tEdT 



< t E >- 



2.1. Pixel lensing basics 



(14) 



Pixel lensing technique is based on the observation of 
the flux variations of every element (pixel) of an image 
l|Ansari et allll997|l . Looking towards M31 a large num- 
ber of stars contribute at the same time to the flux received 
by each pixel so that only highly magnified events can be 
detected. 

To be detectable a microlensing event must give rise to 
a substantial flux variation with respect to the background 
A^bi = -/Vgaj + Ngky, which is the sum of the M31 surface 
brightness and the sky contribution. The excess photon 
count per pixel due to an ongoing microlensing event is 



AA^pix = N hl [A pix - 1] = f see N s [A(t) - 1] 



(15) 



where N s is the source photon count in the absence of lens- 
ing, A(t) is the source magnification factor due to lensing 
see e.g. and /see the fraction of the seeing 

disk contained in a pixel. Therefore, the expected number 
of photons in a pixel will be N p j X = A^i + AiV p j x . 

Of course, a pixel lensing event is detectable if the 
excess pixel photon count is greater than the threshold 
pixel noise ctt- Accordingly, by requiring the signal to be 
at least 3ctt level above the baseline count, one ob t ains a 
threshold value for the amplification ( Kerins et all l200l[) 

}>,n ' (16) 



A 



1 



fsee-Ns 

which corresponds to a threshold value ut for the impact 
parameter, via the well known relation between the lens 
impact parameter and the amplification factor. 

One can estimate or as the maximum between the 
statistical error oc y/TTu and ~ 3 x 10 -3 A^i that is deter- 
mined by the pixel flux stability. Accordingly, ut depends 
on both the line of sight to M31 and the source magnitude 
M. 

Hence, by averaging on the source luminosity function 
4>(M), we can evaluate the average threshold impact pa- 
rameter for any direction towards the M31 galaxy, so that 
we get 

/ u T {x iy ;M) 4>{M)dM 



(u T (x,y))0 = 



(17) 



/ 4>{M)dM 

where the coordinates x and y span the sky plane towards 
M31. 

By using the threshold impact parameter defined in 
ea. Ill 711, one obtains the pixel lensing rate a s follows 



llKerins et alll200lil20o3lDe Paolis eTall l2005h 

F P (x,y) = {urix.y))^ T(x,y) 



(18) 
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MACHO — to — self lensing probability ratio 




Fig. 2. MACHO-to-self lensing probability ratio 
(i\/P s )An map projected onto the sky plane. Here 
and in the following figures we assume a fraction 
/h = 20% of dark matter and a MACHO mass ph. = 0.5, 
both in M31 and MW halos. From the inner to the outer 
part (Ph/f s )An increases with lines referring to the values 
0.5, 1, 2, 3, 4 and 5, respectively. 



MACHO — to — self lensinq Einstein time ratio 




M31 surface brightness and rotation curve profiles, can be 
considered as an acceptable model for the mass distribu- 
tion in the M31 galaxy. Accordingly, the mass density of 
the M31 disk stars is described by a sech-squared profile 



p D (R,z) = p D (0)exp{-R/h) sech 2 (z/#), 



(19) 



where H = 0.3 kpc, h = 6.4 kpc and p D {0) = 0.35 x 
10 9 Mq pc -3 are, respectively, the scale height and scale 
lengths of the disk and the disk central mass density. R 
is the distance on the M31 disk plane (described by the 
coordinates x and y) and z is the distance from it. The 
M31 disk is assumed to be inclined by an angle i = 77° 
and the azimuthal angle relative to the near minor axis is 
cj) = -38.6°. 

The M31 bulge is parameterized by a flattened power 
law of the form 



Pb{R, z) = pb{0) 



a 



(!)' 



-8/2 



(20) 



where pb(0) — 4.5 x 10 9 M kpc -3 , q ~ 0.6 is the ratio of 
the minor to major axis, a ~ 1 kpc and s ~ 3.8. Both the 
M31 disk and bulge are truncated at a distance R = 40 
kpc. 

We remark that the twisting of the optical isophotes 
in the inner M31 regions indicates that the bulge ma- 
j or axis is offset b y ~ 15° from the disk major axis 
l|Stark and Binnevl |l99j). The consideration of this ef- 
fect by Kerins et al. 120051 leads to the evaluation of pixel 
lensing rates that show spatial distributions tilted of the 
same amount inside 5 arcmin from the M31 center. The 
twisting effect vanishes at larger distances due to the in- 
creasing contribution to microlensing by M31 disk and 
halo. Clearly, our results in Figs. 11131 do not show the 
above mentioned effect since here we are considering only 
a flattened bulge without twist. However, we expect that 
the consideration of the isophote twisting does not sub- 
stantially modify our results about both MACHO-to-self 
lensing probability and event time scale ratios (given in 
Tabs. EEJj particularly for events at large distance from 
the M31 center. 

The dark matter in the M31 halo is assumed to follow 
an isothermal profile 



Fig. 3. The MACHO-to-self lensing Einstein time ratio 
< tE h > / < tE s > map is given for p^ = 0.5. Going 
from the inner to the outer part the ratio decreases with 
lines corresponding to values in between 2 to 1.25 with a 
step of 0.25. 



3. Modeling 

3.1. M31 and Galaxy mass distribution models 

The M31 disk, bulge and halo mass distributions are de- 
scribe d adopting the reference model discussed in iKerins! 
( 200'3). This model, providing remarkable good fits to the 



Pff(r) = ph(0)- 



(21) 



with core radius a = 4 kpc and central dark matter density 
p H (0) = 6.5 x 10 7 M Q kpc -3 . The M31 halo is truncated 
at 100 kpc with asymptotic rotational velocity v ro t — 235 
km s . 

We do not consider other dark matter di s tribut ion 
model s, as King-Michie jBinnev and Tremainel Il987|) or 
NFW ilNavarro et allll997^ suggested by N body simula- 
tions. The effect of using such models, which have a more 
concentred dark mass distribution, is both to decrease the 
spatial distribution of MACHOs at large distance from 
the M31 center, where the rotation curve is poorly deter- 
mined, and to increase it in the innermost region, where 
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the MACHO contribution is relatively unimportant with 
respect to that of bulge and disk. Regarding the former 
aspect, we also note that in the MEGA experiment the 
typical MACHO lens distance (about 20 kpc) is too small 
to appreciate the effect of this choice. Overall, the cur- 
rent data do not allow one to perform a meaningfully fine 
tuning of the dark matter parameters. 

As usual, the mass density profile for the MW disk is 
described with a double exponential profile 

Pd(R,z) =p D (Ro)exp(-(R-Ro)/h) exp(-\z\/H) , (22) 

with Earth position from the Galactic center at R ~ 8.5 
kpc, scale height H ~ 0.3 kpc, scale length h ~ 3.5 kpc 
and local mass density pn(Ro) — 1-67 x 10 s M & kpc -3 . 

The dark halo in our Galaxy is also assumed to follow 
an isothermal profile with core radius a ~ 5.6 kpc and 
local dark matter density ph(Ro) — 1-09 x 10 7 M Q kpc -3 . 
The corresponding asymptotic rotational velocity is v rot ~ 
220 km s -1 . The MW halo is truncated at R ~ 100 kpc. 

For both M31 and MW halos, the fraction of dark 
matter in form of MACHOs is assumed to be /h — 0.2 
l|Alcock et al! f20001. However, most of our results can eas- 
ily be rescaled to get the corresponding figures for other 
values of /macho- 

Moreover, we assume that the random velocities of 
stars and MACHOs follow Maxwellian distributions with 
one-dimensional velocity dispersion a = 30, 100, 166 km 
s -1 and 30, 156 km s -1 for the M31 disk, bulge, halo and 
MW disk and halo, respectively. In addition, a M31 bulge 
rotatio nal velocity of 30 km s -1 ha s been taken into ac- 
count l|Kerins et alll200lHAn et all l2004|) . 

3.2. Stellar luminosity function 

Pixel lensing event detection by the MEGA collaboration 
is performed in the red band 1 and, thus, red giants are 
the most luminous stars in this band. Therefore, we may 
safely assume that the overwhelming majority of the pixel 
lensing event sources are red giants. 

Moreover, in the lack of precise information about the 
stellar luminosity function in M31, we adopt the lumi- 
nosity function derived from the stars in the Galaxy and 
assume that it also holds for M31. 

Accordingly, following i Mamon and Soneiralll982|) we 
assume that the stellar luminosity function does not de- 
pend on the position and in the magnitude range — 6 < 
M < 16 it is proportional to the expression 



<p*(M) oc 



1 + io-("-P)S(M-M*)n/s 



(23) 



where, in the red band, M* = 1.4, a ~ 0.74, = 0.045 
and S = 1/3. Moreover, the fraction of red giants (over the 



1 Indeed, observations in the red band, by minimizing light 
absorption in M31 and MW disks by the intervening dust, offer 
the best compromise between sampling and sky background. 
Observations in other bands (B and V) are commonly used to 
test achromaticity of the candidate events. 



total star number) as a func tion of M is approximated as 
llMamon and Soneiralll982ft 

f RG (M) = 1 - Cexp[a(M + (3) 1 ] for - 6 < M < 3 

= for M > 3 , (24) 

where, in the red band, C ~ 0.31, a ~ 6.5 x 10 -4 , j3 — 7.5 
and 7 ~ 3.2. 

Therefore, the red giant luminosity function will be 
</>rg (M) oc (M) x /no (M) and the fraction of red giants 
averaged on the magnitude is given by 



, _ Sl^KG{M)dM 



5.3 x 10 



-3 



(25) 



from which it follows that the local number density of red 
giants will be urg — 5.3 x 10 -3 n*, where n* is the local 
stellar number density. 

3.3. Mass functions 

As concerns the lens mass function ipo(p) m eqs. JHJ) - (|11|) . 
for lenses belonging to the bulge and disk star populations, 
the lens mass is as sumed to follow a broken power law 
llGould et allll997h 



K p- 56 
K p- 2 20 



for 0.1 < (j, < 0.59 
for 0.59 < p < Pup 



(26) 



where the lower limit p m i Yl = 0.1 and the upper limit p up 
is 1 for M31 bulge stars and 1 for M31 and MW disk stars 
(see also Kerins et al., 2001). K is fixed according to the 
normalization given by eq. ©■ The resulting mean mass 
for lenses in the bulges and disks are (mi) ^0.37 M Q and 
{fnd) ~ 0.69 Mq, respectively. 

For the lens mass in the M31 and MW halos we assume 
the (^-function approximation 



^o(m) 



Pa 



M Q p h 



8{p - p h ) 



(27) 



and we take a MACHO mass, in solar units, p^ = 
10- 2 , 10-\ 0.5, 1. 

4. Analytical results 

In the present analysis we adopt the parameters for the 
INT and the Sloan-r filter on the WF C (Wide-Field 
Came ra) used by the MEGA collaboration l)de Jong et all 
120041) . The Telescope diameter, the pixel field of view 
and the image exposition time are 2.5 m, 0.33 arcsec and 
texp = 760 s, respectively. We also use a gain or con- 
version factor of 2.8 e~/ADU, and a loss factor ~ 3, for 
both atmospheric and instrumental. The zero-point with 
the Sloan-r WFC turn s out to be ~ 24.3 mag arcsec -2 
l|Belokurov et al.|.l2005|) . Moreover, we adopt a value ~ 1.5 
arcsec for the average seeing conditions, a sky background 
m s k y — 20.9 mag arcsec -2 (corresponding to a Moon 
eclipse) and a minimum noise level of ~ 2.5 x 10 _3 iVbi. 
iVbi is the baseline photon count which is the sum of the 
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M31 surface brightness given bv lKentl l|l989|) and the sky 
contribution. 

Maps of optical depth, expected event number and 
event timescale in pixel lensi ng experiments have b een 
presented in a previous paper l|De Paolis et all 12005^ to- 
gether with the study of the dependence on microlensing 
quantities with the as sumed M31 mass distribution model 
(see also Kerins EnoiL Our results are also in good agree- 
ment with previous analytical estimates for the rate, the 
timescale distribution and the optical depth l|Baltz et all 
l2003UGvuk and Crottsl Effloh . 

In Fig. ^the map of (ut.(x, y))<f, RG shows that in the 
central M31 regions the lens impact parameter on average 
is < 0.04 implying that only high magnified events with 
At > 25 are in principle detectable. The asymmetrical 
shape is due to the internal extinction of the M31 disk fo r 
which we use the value 0.74 mag given bv iKentJ lfl989). 
Indeed, due to the inclination of the M31 disk, along a line 
of sight towards the southern region there exists a larger 
number of source stars (with respect to the corresponding 
northern field) which are not absorbed by the M31 disk 
dust (and therefore appear with a smaller magnitude), 
lying at larger averaged values for ut- 




M 



Fig. 4. Plot of the source absolute magnitude distribution 
for the generated events (solid line) and for the revealed 
events (dashed line) towards the MEGA 7 direction. 

In Fig. [3 assuming a MACHO mass /ih = 0.5 and 
a halo MACHO fraction /h = 0.2, we show the map of 
the MACHO-to-self lensing probability ratio (Ph/-fs)An- 
We find that in the M31 central regions microlensing is 
dominated by self-lensing contributions while MACHO- 
lensing becomes important at distances > 10 arcmin from 
the center. In the figure it is also evident the well known 



Table 2. MACHO-to-self lensing probability ratio 
(Ph/-Ps)An towards the 14 MEGA events are given for dif- 
ferent MACHO mass values. 



MEGA 


= 0.01 


0.1 


0.5 


1 


1 


1.65 


0.52 


0.23 


0.16 


2 


1.76 


0.56 


0.25 


0.18 


3 


2.22 


0.71 


0.31 


0.22 


4 


2.52 


0.78 


0.35 


0.24 


5 


24.58 


7.72 


3.46 


2.40 


6 


23.46 


7.53 


3.38 


2.40 


7 


12.93 


3.83 


1.80 


1.26 


8 


9.79 


3.15 


1.40 


0.96 


9 


19.49 


6.15 


2.79 


1.99 


10 


17.19 


5.44 


2.42 


1.73 


11 


19.22 


6.04 


2.74 


1.88 


12 


18.91 


5.88 


2.62 


1.91 


13 


54.13 


17.60 


7.87 


5.67 


14 


112.32 


35.47 


17.07 


11.12 



Table 3. MACHO-to-self lensing Einstein time ratio 
(is h)/(tE s) towards the 14 MEGA events are given for 
different MACHO mass values. 



MEGA 


fx h = 0.01 


0.1 


0.5 


1 


1 


0.30 


0.96 


2.15 


3.02 


2 


0.31 


1.00 


2.42 


3.24 


3 


0.29 


0.94 


2.02 


2.91 


4 


0.20 


0.65 


1.47 


2.12 


5 


0.16 


0.51 


1.10 


1.59 


6 


0.16 


0.49 


1.09 


1.55 


7 


0.15 


0.52 


1.10 


1.58 


8 


0.13 


0.41 


0.92 


1.35 


9 


0.13 


0.42 


0.97 


1.34 


10 


0.17 


0.54 


1.21 


1.69 


11 


0.16 


0.50 


1.11 


1.61 


12 


0.14 


0.44 


0.98 


1.36 


13 


0.16 


0.49 


1.08 


1.51 


14 


0.15 


0.48 


1.03 


1.52 



near- far disk asymmetry due to the inclination of the M31 
disk llCrottsl llQflj iBaillon et all Il993t Ijeteerl Il994h . 

In Tab. we show how the MACHO-to-self lensing 
probability ratio (Ph/Ps)An depends on the MACHO mass 
/Lth towards the 14 MEGA events. The general trend is that 
(Ph/-Ps)An increases as /ih decreases as a consequence of 
the increase of the MACHO number density. Moreover, for 
a given value of /ih, (Ph/Ps)An increases with the distance 
from the M31 center. 

In Fig. the map of the MACHO-to-self lensing 
Einstein time ratio (tE hA_E s)An is given for a MACHO 
mass of 0.5 Mq. In the region inside ~ 5 arcmin from 
the M31 center microlensing events due to MACHOs have 
twice as long a duration as compared to self-lensing events, 
while in the regions far away from the M31 center all 
events have roughly the same duration. Obviously, the 
MACHO-to-self lensing Einstein time ratio depends on the 
assumed MACHO mass. This fact is clearly seen in Tab. 
where the ratio (< is h > / < tE s >)An towards the 
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14 MEGA events is given for different MACHO mass val- 
ues, /ih = 0.5 corresponds to the MACHO mass value for 
which the 10 outer MEGA events (5-14) are characterized 
by having a MACHO-to-self lensing Einstein time ratio 
equal to unity, i.e., events due to either MACHO-lensing 
or self-lensing have roughly the same duration. Thus, these 
different types of events are indistinguishable on the ba- 
sis of their timescale alone. The situation is instead much 
more favorable for MACHO masses smaller or larger than 
0.5 M . 

We would like to emphasize that in pixel lensing exper- 
iments ts is not a directly observable quantity, since the 
relevant time scale is the full-width half-maximum event 
duratio n ti /?, which de pends on f B and the impact param- 
eter u llGondoldllflfl^ . However, since the probability for 
a given uq value is practically the same for self-lensing and 
MACHO-lensing events (as we have verified by using the 
MC code) the ratio (< is h > / < tg s >)aii is equivalent 
to the ratio (< t\/i h > / < t\/2 a >)An, at least for the 
MC generated events. Of course, as it will be more clear 
in the following Sections (in particular from Tab. 01 and 
Fig. ll(J|l . it does not mean that the same conclusion holds 
for the MC revealed events, since the (normalized) event 
number which pass the MEGA selection criteria turns out 
to depend on t\/2 and i? m ax (see Fig.EJ). 

5. Monte Carlo simulation 

Once the event location (one of the 14 MEGA events to- 
wards M31) has been selected, for any lens and source pop- 
ulation present along the line of sight we have to make as 
next the choice over the following five parameters: source 
distance D os , lens distance D Q i, lense effective transverse 
velocity v\±, lens mass fj, and source magnitude M. We 
shall denote these parameters by Xi, with i — 1,..,5 in 
the order just listed (e.g. x\ — D os , etc.). The probability 
with which we select the events according to one of the 
parameters is then given by 

P(xi)dXi = -L^^da* (28) 



N ev dXi 



where is defined as being the integrand of the event 
number N ev integrated over all variables Xj without, how- 
ever, the considered variable .r,-. Clearly when integrating 

8N, 



dx over dxi we obtain again N cv . 

This way, for instance, for self-lensing events the prob- 
ability of extracting a lens with mass fi turns out to be 

P{x^ = p)dp cx /i 1 / 2 ipo(fi)dfi , (29) 

where V'o(m) is the lens number density distribution de- 
fined in cq. 

The probability of a source 2 at distance D os and a lens 
at distance D i from the observer, respectively, is given by 




bulge-disk 



Fig. 5. The lens transverse velocity v\± distribution for 
the MC generated events towards the MEGA 1 direction 
is shown in the case of disk-disk (solid line) and bulge-disk 
events (dashed line). 



MC event number density 




Fig. 6. MC revealed event number density Af c ^ cv plot 
(normalized to the maximum value) as a function of 
ti/2 and Rmax- Here we consider both self-lensing and 
MACHO-lensing events, averaged over all the 14 MEGA 
directions and we assume /ih = 0.5 and /h = 0.2. From the 
inner to the outer part Af c ^ cv decreases with lines referring 
to values 0.7, 0.5, 0.3, 0.2, 0.1, 0.05. 



and 



rD os 



P(D os )dD os cx p s (D os )D 3 J s 2 

2 Here we take into account that the source number inside 
the pixel solid angle increases with the distance as D os 2 . 



P(D ol )dD ol dD os (30) P(D ol )dD ol cx Pl [D ol (D os - D ol )Y' 2 (D ol ) dD ol 



(31) 



where p\(D a i) and p s {D os ) are, respectively, the lens and 
the source mass densities. 
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Fig. 7. The i? max distribution for the MC revealed events 
towards the MEGA 14 (solid line), 7 (dashed line) and 1 
(dotted line) directions. 



The probability of extracting a source with magnitude 
M is weighted by the function 



P(M)dM oc u T (M)(f> RG (M)dM 



(32) 



where </>rg(M) is the source magnitude distribution for 
red giants defined by means of eqs. flty - l|2"4*j) and, through 
ut(M), we take into account the surface brightness vari- 
ation along the field of view. Indeed, we find that the me- 
dian value of the source magnitude distribution decreases 
from M median = -2.18 for the innermost MEGA direc- 
tions to 

M median = —1.18 for the outermost ones. In Fig. 
01 we show the distribution of the source absolute magni- 
tude for all MC generated events and for the revealed (see 
after) events towards the direction of the MEGA 7 event. 
For a given source and lens position the probability of ex- 
tracting a lens with transverse velocity v\± is given by eq. 
I)12[l. where z = v\±/<r\. In Fig. for both disk-disk and 
bulge-disk events, we show the lens transverse velocity dis- 
tributions for MC generated events towards the MEGA 1 
direction. The shorter < vi± > value in the first case will 
lead to an increase in the event timescale. 

The probability of extracting an event with impact pa- 
rameter x§ = uq is constant for values of uq in the interval 
in between and ut, as defined in eq. i|17|) 3 - Thus it is 
assumed that 



P(uq) oc const. 



(33) 



3 In order to avoid a sharp cutoff in several obtained plots 
(see, e.g., Fig. [7J we allow for events with «o > «t to be de- 
tected. 



Once all the parameters have been fixed, for each event 
we build the corresponding lightcurve using the same time 
sampling of the MEGA campaign and the same observing 
and instrumental conditions of the considered experiment. 
We add a gaussian nois e modulated by the Moon phase 
and use a Paczyhski fit l|Paczvhskil Il98fj|) to evaluate the 
microlensing parameters. 

The selection of the MC generated events is based 
on the same criteria adopted by the MEGA collabora- 
tion. By using the results of the Paczyhski fit we filter 
the lightcurves with the following selection criteria: peak 
sampling, peak significance, peak width, baseline flatness 
an d goodness of fit (fo r more details see Appendix A.l 
in Ide Jong et alJ 12004^ . The criteria that more severely 
cut the MC generated events are the peak sampling and 
the peak significance, which depends through the evalua- 
tion of the statistical errors on the the Moon phase. The 
former criterium leads to select events with well sampled 
lightcurves (within the MEGA campaign), the latter to 
events with high signal-to- noise ratios. The fraction of 
MC selected events also depends on the event location 
through the threshold value of i? max . This effect is further 
discussed in the following section. 



6. Monte Carlo results 

MC results allow to estimate the features of the re- 
vealed events which pass the adopted selection criteria. 
In Fig. HJ1 assuming a MACHO mass /ih = 0.5 and a halo 
MACHO fraction /h = 0.2 we give the contour plot, in the 
{ti/2, -Rmax) parameter space, of the event number density 



TV rcv = 

•* v PV 



d 2 N. 



dti/ 2 dR„ 



(34) 



averaged on all the 14 MEGA directions, for both self- 
lensing and MACHO-lcnsing events. We also give in the 
same parameter space the position of the 14 observed 
MEGA events. Af ^° v is maximum in the region t\i2 — 20 
day and i? max — 22. Moreover, J\f c * cv rapidly decreases 
for i? max > 23.5 (due to the lack of revealed events with 
low signal-to-noise ratio) and ti/ 2 < 5 day (due to the 
peak sampling) and in the region of high amplification and 
long duration events (upper-left region of the figure) , due 
to the absence of MC generated events. Actually, the last 
cutoff depends on the adopted MACHO mass value and 
shifts towards smaller ti/2 values decreasing Indeed, 
in Fig. Othe crucial parameter determining the event dis- 
tribution is the lens mass value and it turns out that the 
region where Af C v CV is maximum scales with /ih in the same 
way as i x , 2 median and -R max cdian (see Tab. 01 where we give 
the median values of the distributions of the MC revealed 
events as a function of the lens star mass /i», ti/ 2 and 

^max ) • 

From Tab. 01 for /ih = 0.5 one can infer that the Af ^ ov 
plot does not vary substantially for MACHO-lensing and 
self-lensing events occurring away from the M31 center, 
since the lens stellar mass is of the same order of the 
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Table 4. For the MC revealed events (which pass the MEGA selection criteria), the median values ^ mcdlan (column 
2) of the lens star mass distribution for self-lensing events, i 1/2 mcdlan (columns 3-7) and -R ma ™ odlan (columns 8-12) for 
self-lensing and MACHO-lensing are given towards the 14 MEGA events. 





median 

A 4 * 


i median 
I l/2 


T) median 

"max 




self 


self 




MACHO 








self 




MACHO 








MEGA 






/ih = 0.01 


0.1 


0.5 


1 






H h = 0.01 


0.1 


0.5 


1 


1 


0.48 


15.66 


4.96 


11.22 


19 


51 


22.49 


21 


56 


21.16 


21.45 


21 


65 


21.61 


2 


0.48 


15.33 


5.06 


11.22 


18 


15 


22.04 


21 


62 


21.31 


21.52 


21 


64 


21.70 


3 


0.49 


16.38 


5.17 


11.20 


19 


32 


22.62 


21 


69 


21.41 


21.64 


21 


76 


21.74 


4 


0.62 


18.16 


5.04 


11.43 


18 


81 


23.04 


21 


87 


21.57 


21.85 


21 


93 


21.94 


5 


0.46 


20.62 


4.85 


11.43 


19 


80 


23.59 


22 


22 


21.85 


22.11 


22 


23 


22.22 


6 


0.46 


21.43 


4.77 


11.57 


19 


88 


24.27 


22 


17 


21.85 


22.11 


22 


23 


22.23 


7 


0.50 


20.72 


4.95 


10.97 


19 


12 


22.96 


22 


10 


21.79 


22.02 


22 


14 


22.19 


8 


0.65 


21.47 


5.15 


11.60 


19 


90 


23.65 


22 


13 


21.81 


22.05 


22 


18 


22.21 


9 


0.57 


24.04 


5.20 


12.21 


19 


94 


24.17 


22 


19 


21.85 


22.08 


22 


18 


22.20 


10 


0.45 


20.07 


4.73 


10.79 


18 


98 


23.09 


22 


19 


21.83 


22.06 


22 


20 


22.24 


11 


0.46 


20.37 


4.88 


11.02 


19 


79 


23.66 


22 


14 


21.86 


22.08 


22 


21 


22.22 


12 


0.52 


22.30 


5.00 


11.77 


20 


24 


23.94 


22 


17 


21.84 


22.05 


22 


20 


22.21 


13 


0.44 


22.58 


4.98 


12.24 


21 


14 


25.28 


22 


32 


21.99 


22.20 


22 


34 


22.35 


14 


0.43 


24.15 


5.34 


12.58 


21 


90 


26.53 


22 


34 


22.03 


22.25 


22 


35 


22.38 



Table 5. MACHO-to-self lensing probability ratios 
(Ph/P S )MC towards the 14 observed MEGA events are 
given for different MACHO mass values. Probabilities are 
now calculated from eq. H35|) by considering microlensing 
rates, ii/ 2 and R max distributions for the MC revealed 
events. The results in the table scale with the MACHO 
fraction value as /h/0.2. 



MEGA 


fi h = 0.01 


0.1 


0.5 


1 


1 


3.65 


0.70 


0.18 


0.10 


2 


5.71 


0.79 


0.16 


0.08 


3 


6.70 


1.06 


0.23 


0.14 


4 


0.36 


0.56 


0.33 


0.24 


5 


63.51 


13.41 


3.64 


1.98 


6 


7.56 


6.56 


3.48 


2.35 


7 


0.45 


1.83 


1.46 


3.04 


8 


2.67 


2.88 


1.46 


0.97 


9 


109.47 


14.32 


3.58 


1.74 


10 


0.98 


1.95 


1.98 


1.79 


11 


62.95 


8.29 


2.56 


1.28 


12 


0.04 


0.55 


1.36 


1.79 


13 


20.81 


23.20 


9.44 


4.98 


14 


26.99 


29.88 


18.50 


12.24 



MACHO mass. Moreover, for any value of /ih, the Af e ^ ev 
distribution is also weakly dependent on the selected direc- 
tion towards M31. The situation is completely different for 
/ih 7^ 0.5, and this means that self-lensing and MACHO- 
lensing events lie on different regions in the (fi/2, -Rmax) 
parameter space of the corresponding A/ C v CV plot. 

The dependence of t 1/2 mcdian and i? ma ™ odian on the lens 
mass and the MEGA direction is also clear from Tab. 0] 
For self-lensing events, t 1 y 2 medlan increases from the inner 
to the outer part of the M31 galaxy, and in this region 
t 1 medlan also increases for increasing values of ^ mcdlan ; 
following the dependence of the Einstein radius with the 



lens mass 4 . Moreover, for MACHO-lensing events with a 
fixed /ih value, i 1 ^ 2 mcdlan is almost the same for any MEGA 
direction, while it increases, as expected, with the lens 
mass. 

In Tab. 01 is also evident the decrease, moving towards 
the M31 center, of the i? ma ™ edian value that follows the 
surface brightness variation along the field. In particular, 
we find a shift of about 1 mag going from the inner to 
the outer regions. In Fig. [7| we give the (revealed) event 
distributions as a function of R max towards the MEGA 
1, 7 and 14 directions. This clearly shows how the event 
location affects the fraction of the revealed events. 

The event distributions as a function of tin, for both 
self-lensing (solid line) and MACHO-lensing (dashed line) 
revealed events are shown in Fig. |H1 for /ih = 0.5 and for 
the 1 and 7 MEGA directions towards M31 (representative 
of inner and outer events, respectively). Only in the case 
of the MEGA 1 direction the obtained distributions are 
markedly different for self-lensing and MACHO-lensing, 
since self-lensing events are shorter than MACHO-lensing 
events. For the MEGA 7 direction (and the other outer 
directions) self-lensing and MACHO lensing with the same 
lens mean mass have roughly the same t\/2 distributions. 

Our MC results can be used to estimate the lens nature 
and location of the 14 MEGA candidate events, by weight- 
ing the microlensing rate - giving the analytical estimates 
of the MACHO-to-self lensing probability ratio (Ph/P s )An 
shown in Tab.0- with the revealed event number density 
distribution and taking also into account the observed fea- 
tures of the MEGA events. In Tab. assuming a halo 
MACHO fraction / h = 0.2 and different MACHO mass 



4 For self-lensing events, the lens median mass value changes 
with the MEGA direction, since the disk lens mass is on av- 
erage higher than the bulge lens mass and the disk-to-bulge 
probability ratio depends on the MEGA direction. 
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Fig. 8. The t\i% distributions are shown for MEGA 1 
and 7 events and for /i n = 0.5. In each panel, the solid 
curve refers to self-lensing events and the dashed curve to 
MACHO-lensing events. 

values, we give the MC MACHO-to-self lensing probabil- 
ity ratio 

c 7(t ? hs ) \F rev / n obs\ 
h^l/2 1 M PV h {-ti m ^ ) , . 

(JO) 

MC 



M rev < 

cv 

An "cvs 



where we indicate with A/1, 



reV (^max S ) the 



number of events (either MACHO or self-lensing) with du- 



Fig. 9. Plot of (P h /P s ) A n versus (P h /P s ) M c for all the 14 
MEGA events. In panel a) we consider a MACHO mass 
/ih = 0.01 and we use the symbol plus; in panel b) we take 
/ih = 0.1, 0.5, 1 and use symbols time, circle, asterisk, 
respectively. 



ration and magnitude at maximum, respectively, within 2 
standard deviations around the observed values. Here, we 
remark that each Af cv distribution is normalized to the 
total number of revealed events. A comparison between 
the analytical and MC results given, respectively, in Tabs. 
13 and El is presented in Figs. EK and Eh, where we plot 
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Fig. 10. Plot of {t x /2 h/h/2 s )An versus (i 1/2 h/h/i s )mc, 
for = 0.01,0.1,0.5, 1. The symbols are as in Fig. In 
panel a) we consider MEGA 1-4 and in panel b) MEGA 
5-14 events, respectively. 



(P h /P s )An versus (P^/P^mc, for the MACHO mass val- 
ues fih = 0.01 and /ih = 0.1, 0.5, 1 , respectively. Only 
for /ih = 0.5 and for the outer events MC and analytical 
estimates are in good agreement, since, as we have already 
discussed above, MACHO-lensing and self-lensing events 
have for this particular mass value, on average, the same 
features. 



Also for the time scale ratios between MACHO-lensing 
and self-lensing events analytical and MC results give 
different estimates. A comparison of the analytical re- 
sults in Tab. EH and MC results in Tab. 0] is presented in 
Figs. 1 10b and 1 10b. where we show {ti/2 hAi/2 b)ah ver- 
sus (tx/2 hAi/i s)mc for MEGA 1-4 and 5-14 events. The 
difference is particularly important for the events 1-4. 

It emerges, therefore, clearly that the MC analysis is 
essential for determining the lens nature and location of 
microlensing events, at least if the MACHO mass dif- 
fers substantially from the average self-lensing mass. From 
Tab. 03 we find that MEGA events 1-4 are most likely self- 
lensing events for MACHO masses greater than 0.1 solar 
masses, while events 5, 6, 9, 11, 13, 14 are likely MACHO- 
lensing. For the other MEGA directions the lens nature 
is more uncertain. As a final comment on Fig. we note 
that the event 12, lying in a region of low (revealed) event 
number density, hardly can be considered a reliable mi- 
crolensing event, unless the MACHO mass is considerably 
greater than 1 M Q . 

7. Conclusions 

We have studied the main features of the expected mi- 
crolensing events in pixel lensing observations towards 
M31, by using both analytical estimates (from the mi- 
crolensing rate) as well as results by a MC code where 
we reproduce the observing and instrumental conditions 
of the MEGA experiment. 

First of all, we derive in Section 2 the microlensing 
rate, and assuming a specific mass distribution model 
for M31 and the Galaxy, we calculate the MACHO-to- 
self lensing probability and the MACHO-to-self lensing 
event time scale ratios. For /ih > 0.1, we find that self- 
lensing dominates in the M31 central regions. Moreover, 
for /Xh — 0.5, towards the innermost MEGA directions, 
MACHOs events have duration twice as long as self- 
lensing events, while outer events have roughly the same 
duration. 

We then generate a large number of MC microlensing 
events by choosing relevant source and lens parameters 
as outlined in Section 5. We study the observability of 
the MC events, by considering the capabilities of the INT 
Telescope, typical observing conditions and the event se- 
lection criteria adopted by the MEGA Collaboration. 

MC results can be used to evaluate, for the 14 MEGA 
candidate events, the MACHO-to-self lensing probabil- 
ity and the event time scale ratios (Ph/P S )~MC and 
(ti/2 hAi/2 s )mc (given in Tabs. 4 and 5), by taking into 
account not only the analytical expectations from the mi- 
crolensing rate (as already done in Tabs. 2 and 3) but 
also the features of the MC (revealed) events and the ob- 
served values of ti/ 2 and i? m ax- MC results and analytical 
expectations are compared in Figs. 9 and 10, where one 
can see that in determining the lens nature and location 
of the MEGA candidate events, the MC analysis is partic- 
ularly important for /ih ^ 0.5. Accordingly, we find that 
event 12, lying in a region of low event number density, 
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hardly can be considered a reliable microlensing event, 
unless the MACHO mass is considerably larger than 1 
Mq. Moreover, for a MACHO mass greater than 0.1 M , 
the innermost MEGA events 1, 2, 3, 4 are most likely self- 
lensing events, while 5, 6, 9, 11, 13, 14 are MACHO-lensing 
events. For the other MEGA directions the lens nature is 
more uncertain. 
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